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Abstract 

A nonlinear Schwinger-Dyson (SD) equation for the gauge boson propagator 
of massless QED in one time and two spatial dimensions is studied. It is shown 
that the nonperturbative solution leads to a nontrivial renormalization-group 
infrared fixed point quantitatively close to the one found in the leading order 
of the 1 /N expansion, with N the number of fermion flavors. In the gauged 
Nambu-Jona-Lasinio (GNJL) model an equation for the Yukawa vertex is 
solved in an approximation given by the one-photon exchange and an ana- 
lytic expression is derived for the propagator of the scalar fermion-antifermion 
composites. Subsequently, the mass and width of the scalar composites near 
the phase transition line are calculated as functions of the four-fermion cou- 
pling g and flavor number N. The possible relevance of these results for de- 
scribing particle-hole excitations, in particular antiferromagnetic correlations, 
observed in the underdoped cuprates, is briefly discussed. 
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I. INTRODUCTION 



Quantum electrodynamics in 2 + 1 dimensions (QED 3 ) has attracted much interest 
over recent years. Its version with N flavors of massless four-component Dirac fermions 
shares a number of features, such as confinement and chiral symmetry breaking, with four- 
dimensional quantum chromodynamics (QCD). The loop expansion of a massless theory 
suffers from severe infrared divergencies. However, in the 1/N expansion, the theory be- 
comes infrared finite M, with the effective dimensionless coupling 

e 2 e 2 N /—? 

" (p) = P (i + iW n(p) = ^' p = Vp ' (1) 

giving rise to the renormalization-group (3 function: 

&(a)=p^ = -a(l-^a). (2) 

In Eq. (JTJ) e is the dimensionful gauge coupling and II (p) is the polarization operator. At large 
momenta (p ^> a = e 2 N/8) the effective coupling ([]]) approaches zero (asymptotic freedom) 
while for small momenta (p <C a) it runs to the infrared (IR) fixed point 8/N. Here, the 
dimensionful parameter a plays a role similar to the Aq C d scale. Since QED 3 is a super- 
renormalizable theory, the running of the coupling should be understood as a Wilsonian 
rather than Gell-Mann-Low type, and it is not associated with ultraviolet divergencies. 

By studying the Schwinger-Dyson (SD) equation for the fermion self-energy in leading 
order of the 1/N expansion, it was found in Ref. that a phase transition occurs when the 
coupling at the IR fixed point exceeds some critical value (8/N > 7r 2 /4). This means there 
exists a critical number of fermions N CI (N CT = 32/ it 2 ~ 3.24) below which dynamical mass 
generation takes place and above which the fermions remain massless. This is similar to what 
happens in quenched QED 4 HI], where the gauge coupling must exceed a critical value for 
chiral symmetry breaking to occur. The appearance of a dimensionless critical coupling can 
be traced to the scale invariant behavior of both theories. The scale invariance of QED 3 is 
associated with the IR fixed point, since, as is evident from Eq. ([]]), in the limit p<Ca the 
dimensional parameter e drops out of the running coupling (as well as from SD equations for 
the Green's functions). Related to this is the fact that the chiral symmetry breaking phase 
transition in both theories belongs to a special universality class called conformal phase 
transition (CPT) introduced in Ref. ||. It is characterized by a scaling function having an 
essential singularity at the transition point, and by abrupt change of the spectrum of light 
excitations as the critical point is crossed (for details about the CPT in QED 3 see Ref. f|). 

The presence of a critical N CT in QED 3 is intriguing especially because of possible ex- 
istence of an analogous critical fermion number Nf = N CI in (3 + l)-dimensional SU(A' r c ) 
gauge theories, as is suggested by both analytical studies and lattice computer simu- 

lations |p],[T0[| . Also, a nontrivial IR fixed point in QED 3 may be related to nonperturbative 
dynamics in condensed matter, in particular, dynamics of non- Fermi liquid behavior flI^ , |T2"f . 

The fact that the value of the IR fixed point determines the critical N CI , below which 
the system is in the symmetry broken phase, and that this critical value is found to be 
of order 3 provides motivation for searches beyond the 1/N expansion. It is especially 
important because there is still controversy concerning the existence of finite N cr in QED 3 ; 
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some authors argue that the generation of a fermion mass occurs at all values of N |i3|J14}1 
what might mean the absence of the IR fixed point for the running coupling.]]] Despite 
the arguments of Ref. [15] in addition to the fact that studies of 1/N 2 corrections to the 
gap equation showed the increase of the critical value (N CT = 128/37T 2 ~ 4.32) Hl6|-[T8|], the 



situation is far from being conclusive. What we need is some kind of self-consistent equation 
for the running coupling which is to be solved nonperturbatively. 

In the present paper we study such a nonlinear equation for the running coupling which 
is the analogue of the ladder approximation for the fermion propagator.^] Similar to the 
gap equation, the kernel is taken in the 1/N approximation, where it is nothing else as the 
one-loop photon-photon scattering amplitude with zero momentum transfer. The equation 
obtained is obviously gauge invariant. We then study our equation both analytically and 
numerically. We find that the vacuum polarization operator, obtained as a nonperturba- 
tive solution of the equation, has the same infrared asymptotics as the one-loop expression: 
n(p) ~ Ca/p, C ~ 1 + 1/1AN . Thus a nontrivial IR fixed point persists in the nonperturba- 
tive solution. Moreover, the correction to the one-loop result (C = 1) is small even at iV = 1 
due to smallness of the numerical coefficient before 1/N, that explains why the leading order 
in the 1/N expansion (the one- loop approximation) for the vacuum polarization works so 
well. 

Further we proceed to studying QED 3 with additional four-fermion interactions (the 
gauged Nambu-Jona-Lasinio (GNJL) model). Such kind of models are considered to be 
effective theories at long distances in planar condensed matter physics, in particular, for 
high temperature superconductivity J20|. It is well known that in the improved ladder 



approximation (with the photon propagator including fermion one-loop effects) this model 
has a nontrivial phase structure [f2l| in the coupling constant plane (1/N,g), where g 




2GA/-K 2 is the dimensionless four-fermion coupling (A is the ultraviolet cutoff). The critical 
line is 

1 / / AT \ 

N > N cr , (3) 

at g > 1/4, and 1/N = 1/N cr at g < 1/4. Above this line the gap equation for the 
fermion self-energy has a nontrivial solution. Thus the chiral symmetry is dynamically 
broken, which implies the existence of a nonzero vacuum condensate (ipip)- One end point 
(1/N — 0,g — 1) of the critical line corresponds to the ordinary NJL model (in 2 + 1 
dimensions ), while the other one (1/N — l/N CT ,g = 0) corresponds to pure QED 3 . 

A nice feature of this model is that it is renormalizable in the 1/N expansion leading 
to an interacting continuum (A — * oo) theory near a critical scaling region (critical curve) 



This would happen, for example, if one finds more soft behavior of the polarization operator in 
the infrared, like II(p) ~ (a/p) 7 , with 7 < 1. 

2 Recently, in Ref. Jn| another nonlinear equation for the running coupling was proposed in order to 
study nontrivial infrared structure of the theory. However, their definition of the running coupling 
deviates considerably from the standard one used in present paper and we will not attempt to 
compare both approaches. 
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separating a chiral symmetric phase (xS) and a spontaneous chiral symmetry broken phase 
(SxSB). The spectrum of such a theory contains pseudoscalar (n) and scalar (cr) bound 
states which become light in the vicinity of the critical line. Since the phase transition is 
second order along the g > 1/4 part Eq. (|3]) of the critical curve, scalar and pseudoscalar 
resonances are to be produced on the symmetric side of the curve, whose masses approach 
zero as the critical curve is approached ||23l , |24|j . The part of the critical curve with g < 
1/4, 1/N — 1/N CT is rather special and is related to the CPT in pure QED 3 (we shall discuss 
it more in the main text). 

In this work we study scalar composites (a and tc bosons) which are resonances in the 
symmetric phase of the 2 + 1-dimensional GNJL model. The 7r boson can be viewed as a 
Goldstone boson precursor mode that comes down in energy as the transition is approached. 
Our study is motivated partially by possible relation of these resonances to spin excitations 
observed in neutron scattering experiments in underdoped high-T c superconductors |25 



We calculate their masses and widths as a function of the four-fermion coupling g and 
therefore mass and width's dependence on the doping concentration (since in certain low- 
energy effective models based on spin-charge separation, the coupling g would depend on 
the doping, e.g., Ref. [p6fl ). 

The plan of the present paper is as follows. In Sec. [TT] we derive a nonlinear equation 
for the effective running coupling in pure QED 3 which is then solved both analytically 
and numerically to establish the existence of a nontrivial IR fixed point. In Sec. |T| after 
introducing the GNJL model in 2 + 1 dimensions we solve the equation for the Yukawa 
vertex with nonzero boson momentum. In Sec. IV we obtain an analytical expression for 
the boson propagator valid along the entire critical line and analyze its behavior in different 
asymptotical regimes. The analysis of the scalar composites near the critical line (^|) is given 
in Sec. V. We present our summary in Sec. [VI|. In Appendix [A] we compute the one-loop 
photon-photon scattering amplitude with zero transferred momentum and list some useful 
angular integrals. In Appendix [FJ, an expression for the nonlocal gauge for the ladder (bare 
vertex) approximation is derived. Finally, in Appendix 0, we present some details of the 



approximation [24] which is used to solve the equation for the Yukawa vertex. 



II. THE EQUATION FOR THE RUNNING COUPLING IN QED 3 

The Lagrangian density of massless QED 3 in a general covariant gauge is given by 

C = ~F% + Y a ^ A " )2 + $<*f D id>*> ( 4 ) 

where = <9 M — ieA^ is the covariant derivative. In a parity invariant formulation we 
consider N flavors of fermions (i = 1, . . . N) described by four-component spinors. The 
three 4x4 7-matrices are taken to be 

v=(?_y. ^(Ti,). < s > 

with Oi the usual Pauli matrices. There are two matrices 
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that anticommute with 7 , 7 1 and 7 2 . Therefore for each four- component spinor, there is 
a global U(2) symmetry with generators J, 4-7 3 , 7 s , and |[7 3 ,7 5 ], and the full symmetry 
is then XJ(2N). In what follows we shall restrict ourselves to the symmetric phase of the 
model, i.e., massless fermions. 

The exact SD equations are given in Fig. [I]. For clarity we have extracted the explicit 
factors of iV coming from the one-fermion loop. Since in pure QED 3 we have only one 
dimensionful parameter e, this enables us to choose our scale such that Ne 2 remains fixed. 
This means that every photon propagator (times e 2 ) contributes one factor of 1/N. 

To make a 1/N expansion of Fig. [T], we first need to expand the two-fermion, one-photon 
irreducible fermion-fermion scattering kernel, see Fig. |^. We can convince ourselves that 
Fig. |2| is indeed the right expansion, since the only corrections of order one are fermion 
loops and they are already included in the full photon propagator. Inserting this expansion 
into the SD equation for the vertex, we obtain a closed set of integral equations. The nice 
feature of this truncated system of SD equations is that it satisfies the Ward-Takahashi 
(WT) identities for the vertex as well as for the vacuum polarization . However, finding 



an analytic solution seems to be a formidable task and further approximations are required. 
For such an approximation, we simplify the fermion and the vertex SD equations by keeping 
only the lowest order terms in the 1/N expansion (see Fig. |JJ). Then, inserting this into the 
SD equation for the photon propagator we obtain the equation shown in Fig. |] which still 
satisfies the WT identity. After that equation has been solved, the fermion propagator and 
the vertex can be evaluated explicitly through right-hand sides of Fig. |3|. 

We could solve the obtained photon propagator equation by further iteration, with the 
one-loop fermion correction included at the initial step to obtain a perturbative 1/N expan- 
sion. Instead we choose to solve the nonlinear integral equation given by Fig. |] as it is. In 
this way we might get a hint of any nonanalytic behavior in 1/N which would be lost other- 
wise. At first glance, this way of solving a truncated system of SD equations ignores possible 
nonanalyticity in 1/N coming from the fermion wave function and the vertex (related, for 
example, to the power-law behavior due to an anomalous dimension). Note, however, that 
the fermion propagator is a gauge dependent quantity, thus possible power-law behavior of 
the fermion wave function must cancel the corresponding behavior coming from the longi- 
tudinal part of the vertex (recall that we consider in this paper the symmetric (massless) 
phase only). Therefore the only nonanalyticity we have neglected is the one which might 
be present in the transverse part of the vertex beyond order 1/N. Neglecting possible non- 
analyticity in the transverse vertex means that we are seeking for nonanaliticity originating 
from the nonlinear equation for the photon propagator only. To some extent, the considered 
approximation is similar to solving the SD equation for the fermion mass function in the 
ladder approximation, where the photon propagator is taken in the leading 1/N order 

In 2 + 1 dimensions, the SD equation for the photon propagator reads 

d; u \ p ) = d ^(p) + u^ p ), (7) 

with Dq^v the bare photon propagator, and where is the vacuum polarization tensor 

W v (p) = iNe 2 / — -^3 Tr [j"S{r + p)T u {r + p, r)S{r)) . (8) 
Because of the gauge symmetry the vacuum polarization tensor is transverse: 
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n^(p) = (-^y+ P v)n(p), (9) 

therefore, for the full photon propagator in a general covariant gauge, we have 
n f \ ( , P^A 1 1 PmP" 

= + IT J f ITTnM " %^ (10) 

Moreover, one can write 

n(p) = -^(^-c 1 ^)n^( P ), (n) 

where the constant c± can be chosen arbitrarily. 

The vacuum polarization Tl(p) governs the running of the dimensionless gauge coupling. 
Now we study the integral equation based on Fig. |], this gives 

if» = ny (p) + nr (p) + o (i/ao , (12) 

where II]' 1 ' (p) is the one-loop vacuum polarization, 

nf(p) = ^e 2 / ^Tr[ 7 ^o(r + p) 7 ^ (r)], (13) 

with Sq(p) the bare fermion propagator, So(p) = 1/p, and 

t^^WS^Cp,*), (14) 

where B tJif>ua (p, k) is the one-loop "photon-photon" scattering amplitude, with zero momen- 
tum transfer, i.e., 

B^{p, k) = i — ^- Tr [fS (r + p)i p S {r +p + k)-fS (r + k)-fS (r)] 



m (2tt) 



+ l / 7Tr4Tr[7^o(r + p)7^o(r+p + A;)7 CT So(r + p)7^o(r)] 



+ i ^-T T [^S (r + p)YSo(r)YSo(r + k)YS (r)]. (15) 

jm i^) 

A graphical representation of the "box" diagram ( TTSj ) in terms of Feynman diagrams is given 
in Fig. §. 

For the scattering amplitude B wva there exists a Ward-Takahashi identity WA , which 
states the transversality of the amplitude with respect to external photon momenta, 

B^ pa (p,k,q,r) = 0, k v B^ pa (p, k,q,r) = 0, etc. (16) 

The vacuum polarization tensor has a superficial linearly divergent part, which can be 
removed by a proper gauge-invariant regularization. However, since the divergent part is 
proportional to g pu we can project out the finite vacuum polarization by contracting IT^p) 
with the projector 
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V 



(17) 



i.e., we choose the constant c\ in (O) to be c\ = 3. This approach was used in Refs. |2S 

(18) 



and |15| . In this way, we obtain 

n(p) = n 1 (p) + n 2 (p) + o(i/iV) 



with 



ni(p) 
n 2 (p) 



4iiVe 2 



iNe 4 
~2p^~ 



d 3 k 
(2^)3 



k 2 - 2k-p-3(k-p) 2 /p 2 



m y^r L k 2 (k + p) 2 

d 3 k B(p 2 ,k 2 ,p-k) 

(2^ p [i + n(fc)] ' 



where 



B( P 2 ,k 2 ,p-k)^g fll/ g pa B^{p,k). 
In Euclidean formulation the above expressions can be written as 



n 2 (p) 



2iVe 2 
iVe 4 



r A; 2 - 2k-p-3(k-p) 2 /p 



dQ 
Air 

An 2 pJ ~" p[l + IL(k)] 



.21 



p(k + pY 



Ne 2 
~~8p 



dk 



where 



K{p,k) 



dVL 
An 



B(-p 2 ,-k 2 ,-p- k). 



(19) 
(20) 

(21) 

(22) 
(23) 

(24) 



From Figs. |] and [^, one can see that the first term in Eq. (|i~5|) corresponds to a vertex 
correction and the last two terms are fermion self-energy corrections. The sum of these 
diagrams has symmetries which provide a consistency check on the final result. From the 
graphical representation it is obvious that the quantity B(p 2 ,k 2 ,p ■ k) should be invariant 
under p <r-> k and under p — > —p or p ■ k —>■ —p ■ k. 

A detailed computation of the "box" function B is presented in Appendix [A], and the final 
expression for B is given by Eq. (|A11|) . One can verify that Eq. (|A11|) has the symmetries 
we mentioned above. Finally, we perform the angular integration to obtain K(p, k), 



k ■ p 



1 1 

k p 2kp\k — p 



2k 2 + kp + 2p 2 2k 4 + 5k 2 p 2 + 2p 4 



kp 



1 1 

k p 6max(/c 3 ,p 3 
2k 4 + 5k 2 p 2 + 2p 4 



sinh 



kp\k — p 

2k 2 + kp + 2p 2 

kp max(/c,p) 
_ x min(p, k) 



2{k ■ p)kp\k — p\ 



2k 2 p 2 y/p 2 + k 2 '"' max(p,k) 
where we have made use of the integrals given in Appendix [A]. 



(25) 
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Thus, we arrive at the following nonlinear equation for the vacuum polarization: 



TT/ , Ne 2 Ne 4 f „ K(p, k) , , 

n (p ) = — / dk . v „ /,.. . 26 

w 8 P 47T 2 p7 P [i + n(fc)] 1 ; 

Apparently, this equation is gauge invariant. We can rewrite it also as the equation for the 
running coupling d(p) which must be self-consistently determined from it: 

oo 

N f 

a-\p) = a^(p) -^T J dkkK(p,k)a(k), (27) 

o 

where a,\{p) is the one- loop running coupling (see Eq. ([!])). 

Equation fl27|) is the simplest nonlinear equation for the running coupling (or the photon 
propagator) which is derived at the lowest order in the 1/N truncation of the SD equations. 
In fact, it should be considered as an analogue of the ladder approximation for the fermion 
propagator. The effects of a constant fermion mass can be incorporated at one's wish by 
computing the box diagrams with massive fermions. This would allow one to study the 
coupled system of the SDE for the fermion self-energy and photon polarization operator 
along the lines of Ref. [|18|| . However, this is beyond the scope of the present paper and we 
shall leave aside this issue. f\ 

Now we proceed by solving Eqs. (|2B|) and ( P7| ) both analytically and numerically. Ap- 
proximating, as usual, the expression ( |25| ) for the kernel by its asymptotics at p ^> k and 
p <C k 

K{p,k)^-± P>t (28) 
15 maxyr , k') 



one can reduce the integral Eqs. (|2q ) and (|27| ) to differential ones in order to study the 
asymptotical behavior of H(p) and a(p) in the ultraviolet and infrared regions. However, in 
the present case we can find corresponding asymptotics directly from the integral equations. 

First of all, we can immediately see that the solution of Eq. (|2"T| ) for the running coupling 
possesses a nontrivial IR fixed point. Indeed, by making a change of variables, k — > kp, in 
the integral and assuming that a(0) ^ we come to the quadratic equation for a(Q): 



n r 

cT l {ti) = ar (0) - — ^ / dkkK(l,k)a(0), (29) 
47H J 







where we have made use of the fact that pK(p, kp) = K(l, k), see Eq. (^5|) . The last integral 
can be evaluated exactly (see Appendix 0), and we obtain 



3 A coupled system of SD equations for the vacuum polarization and the fermion renormalization 
wave function was studied in Ref. [15] using an Ansatz for the full vertex satisfying the Ward- 
Takahashi (WT) identity. Though such an approach reproduces standard value for the critical 
./V cr ~ 3.3, it does not permit us to identify the Ansatz with a class of Feynman diagrams. 
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m = m , c, = 2 + 2\/ 1 + ]v^- 2 J= !l + 7v^- 2 J- 1+ i«7v- < 30 » 

This result illustrates that the 1/N expansion is reliable even for a rather low number of 
flavors, e.g. N = 2, because of the smallness of the numerical coefficient in front of the 1/N 
term. 

The next term in the expansion of a(p) at small p can also be calculated exactly, as well 
as its asymptotics at large momenta but we focus on finding the asymptotics of the vacuum 
polarization operator itself. For it we seek a power solution (~ (p/a) 1 ) in both asymptotic 
regions, (p <C a) and (p a). We find that the power exponent can only be 7 = — 1 in 
both cases. Thus we get 

01 

n(p) = C-, for p < a, (31) 
p 

(X 

II(p) = -, for p»a (32) 

with the constant C defined in Eq. ( |3D| ) (we recall that a = e 2 N/8). Hence for the running 
coupling we have 

° M = Ri fca/p) » 4' P<<a ' (33) 
e 2 

a(p) = rr> P > «• ( 34 ) 

p{l + a/p) 



The numerical solution of Eq. (|26| ) is presented in Fig. |6|. From this figure it is clear that the 
IR behavior (i.e., p a) of pll(p) is indeed constant and in agreement with the analytic 
analysis. 

For studying effects like symmetry breaking and dynamical mass generation, it is suffi- 
cient to consider only momenta less than a. Therefore for the remainder of this article we 
will just use Eq. (j3~lD and treat a as the ultraviolet cutoff for nonperturbative dynamics. 
This allows us to write the gauge boson propagator as (in Euclidean formulation) 

e 2 DM = ("^ + (1 - m P ~f) ^f; «(0) = JL (35) 

for p = a/p 2 " < a, with C given by Eq. (0), and where parameterizes a nonlocal gauge 
fixing function (see Appendix fB|). This form of the photon propagator will be used in the 
next section. 

The gauge boson propagator of Eq. (^Sj) gives rise to a Coulomb potential instead of 
a logarithmically confining potential. The dimensionless coupling a = a(0) should now 
be interpreted as the coupling parameter of a perfectly marginal (or conformal invariant) 
interaction: /3(«o) = 0. 
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III. QED 3 PLUS FOUR-FERMION INTERACTIONS 



The gauged NJL model with N fermion flavors is described by the Lagrangian 

G 

4""" ' rv ~ ' ~ M "~"' T ' 2N 



Cgnjl = — A Flu + H^D» ~ ™o)V> + — T [{^f + (36) 



where = — ieA^ is the covariant derivative, and the last term is a chirally invariant 
four- fermion interaction with G the corresponding Fermi coupling constant. In the absence 
of a fermion mass term m which breaks the chiral symmetry explicitly, the Lagrangian 
(|36|) possesses a U(l) gauge symmetry and a global Ux,(1)xUr(1) chiral symmetry. For the 
four-fermion coupling we introduce the dimensionless coupling constant g = 2GA/it 2 , and 
we consider the dimensionful gauge coupling e 2 as the UV cutoff (more precisely, a ~ A). 

A parity invariant bare mass term moipip as well as a dynamically generated fermion mass 
breaks the global symmetry down to XJl+r(1)- Further we study mainly the chiral symmetric 
case with mo = 0. By introducing the auxiliary scalar fields a and 7T, the Lagrangian (36) 
can be rewritten as 

C 2 = ~F, U F^ + ji-fD^ - ^{a + i 7B 7r)^ - ^ (a 2 + rr 2 ) , (37) 

where a = -(G/N)^, vr = -(G/A)^7s^- 

The propagators of the a and 7r fields, As and Ap, are defined, respectively, as follows: 

A s (q) = -ij d 3 x e iqx (0\T(a(x)a(0))\0) G , (38) 

A P (q) = -i J d 3 x e^(0|T(7r(x)7r(0))|0) c , (39) 

where the subscript C stands for "connected." The SD equation for the scalar (pseudoscalar) 
propagator is given by 

a 5(p)(p) = -§ + n 5 (P)(p), (40) 

where the (pseudo) scalar vacuum polarization is 

f A d 3 k 

n s{ P)(p) = t J j^TT[s(k + p)r siP) (k + p,k)S(k)r 0S{P) ] (4i) 

(see Fig. ^), S(k) is the full fermion propagator (5 _1 (A;) = kA(k) — B(k)), and T S (p){k+p, k) 
is the fermion-antifermion (Yukawa) vertex (the bare Yukawa vertices are given by Tqs = 1, 
Top — *75, where 1 is the identity matrix). The absence of kinetic terms for the a and 
ii fields in the Lagrangian is reflected in the constant bare propagator —G. The Yukawa 
vertices Ts and Tp are defined as the "fully amputated" vertices, 



S(k)r s (k,p)S(p)A s (k-p) = -J d 3 xd 3 y e lkx -^(0\T^(x)^(y)a(0))\0) c , (42) 
S(k)T P (k,p)S(p)Ap(k -p) = - [ d 3 xd 3 y e lkx - ip y(0\T^(x)^(y)7r(0))\0) c , (43) 
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In the symmetric phase of the GNJL model the pseudoscalar and scalar propagators are 
degenerate, so are the pseudoscalar vertex and scalar vertex. 

We shall study the SDE for the Yukawa vertex T s and scalar propagator A s with both 
the gauge interaction and the four-fermion interactions treated in the leading order of the 
1/N expansion. This approximation is obtained by replacing the Bethe-Salpeter kernel 
K by planar one photon exchange graph with the photon propagator given by Eq. fl35|) 
and bare fermion-photon vertices (see Fig. |2|). In principle the Bethe-Salpeter kernel also 
contains scalar and pseudoscalar exchanges. One can question whether such exchanges 
can be neglected. In fact, if one includes the ladder like one-scalar and one-pseudoscalar 
exchanges in the truncation of the BS kernel K in the SDE for the Yukawa vertices, then 
such contributions cancel each other exactly in the symmetric phase. On the other hand, 
in the equation for the fermion wave function A(p) these contributions add and must be 
taken into account. Since we take the bare vertex approximation, we need to set A(p) = 1 
for consistency with the WT identity. In Appendix [B| we prove the existence of such a 
nonlocal gauge for the GNJL model in the bare vertex approximation and in arbitrary 
dimensions.^ There it is shown also that four-fermion contributions into the gauge function 
£ are suppressed leading to = 2/3 (Nash's nonlocal gauge). In what follows we use the 
Nash gauge for the photon propagator flSoT). 

The equation for the Yukawa vertex, within the proposed approximation, reads 

r S (p + q,p) = l + ie 2 J — j x S(k + q)T s (k + q, k)S(k)j° D Xu (k - p) (44) 



(see Fig. |J). In the symmetric phase, the equation for the scalar vertex, Eq. (f44|), is a self- 
consistent equation if one uses a gauge where the full fermion propagator has the form of 
the free or bare fermion propagator S(p) = S (p) — 1/p. 

The invariance under parity and charge conjugation restricts the form of the Yukawa 



vertices to the following decomposition [|30| , p4}1 

F S (p + q,p) = l[Fi(p + q,p) + (qp-pq)F 2 (p + q,p)}, (45) 
Fp(p + q,p) = (175) [Fi(p + q,p) + (qp - pq) F 2 (p + q,p)] (46) 

in the symmetric phase. The two scalar functions Fi are symmetric in the fermion momenta: 

F t (p + q,p) = F t ((p + q) 2 ,p 2 , q 2 ) = F t (p 2 , (p + q) 2 , q 2 ), i = 1, 2. (47) 

This is analogous to the four- dimensional case. 

Since we are considering the symmetric phase, the a and tc propagators are identical. 
In what follows we neglect the contribution of F 2 to the Yukawa vertices. The validity of 



this approximation was argued in Ref. [p4J] for the four- dimensional case and the analysis 
can be generalized straightforwardly to the three-dimensional case. Here we only point out 
that calculating F\,F 2 in 1/N perturbation theory reveals that the function F\ contains 



4 A version of a nonlocal gauge in D = 4 leading to approximate equality A = 1 was proposed in 
Ref. p|. 
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logarithmic terms which build up into the power-law form of the full solution (see below); 
on the other hand, F 2 does not contain such logarithmic terms and thus will not contribute 
to the leading and next-to- leading in 1/N order behavior of 

Hence, neglecting all functions except F\, we obtain (in Euclidean formulation) after 
substituting Eq. (B3) with f (p) = 2/3 in Eq. (M) 



A 

F l (p + q ,p) = l + xjdkf§ «•* > I F,(t + «, k), (48) 



where A = 32/(3iVC7r 2 ) and where f dfl denotes the usual angular part of the three- 
dimensional integration. The equation for the a boson vacuum polarization is 

/ x 2iV } „ f dtt (k 2 + q-k) „ 



To resolve the angular dependence of the Yukawa vertex function F\ it is convenient to 
use an expansion in Legendre polynomials P n (see also Appendix 0), 

00 

Fi(p + q,p) = Fi(p,p + g) = f n (p,q)P n (p • g/pg), (50) 

n=0 

where in the right-hand side expression p = a/p 2 , g = a^/^, and p ■ q/pq = cos a. Then we 
follow the arguments of Ref . |24j] and assume that the Yukawa vertex function F\ (p + q, p) 
depends only weakly on the angle p ■ q/pq between fermion and a boson momenta, so that 
the set equations for /„ reduces to the equation for the zeroth-order Legendre coefficient 
function f only This is equivalent to approximating T s by its angular average 

— Ft(p + q,p) = l/ (p,g), (51) 

where the function /o(p, q) depends on the absolute values of the vectors p, q, i.e., in it 
p = \/p*, q = ^/q 2 - Accordingly we write 

/ (p, q) = F m (p, q)9(q - p) + F uv (p, q)9(p - q), (52) 

where the functions Fir and Fuv satisfy integral equations which are given in Appendix y 
(see Eqs. ( |C16|) and (|C17| )). Within this approximation, we find that the scalar vacuum 
polarization (03) is expressed through the function Fuv (see Eqs. (|C9|) , (|U11|), and (|CJ13|j) : 



n S (q) = ^j [F uv (A,g)-l]. (53) 



The integral Eqs. (|C16|) and ( |C17| ) can be reduced to second order differential equations 

tViR + Zp-^Fm + = 0, (54) 



V^F VV + 2p-^-F uv + A ( 1 - £ ) F uv = 0, (55) 



dp 2 dp \ 2p 
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with four boundary conditions. The infrared and ultraviolet boundary conditions (IRBC 
and UVBC), respectively, are 



0. 



p=0 



F\jv(p, q) + P-^F VY (p, q) 



1. 



p=A 



There is a continuity and differentiability equation at p = q: 



Fm{q,q) = F uv (g,g), 



d 
dp 



p=q 



d 
dp 



Fvv{p,q) 



(56) 



(57) 



p=q 



The equation for Fuv can be written as 

-^F UV+ (A-, 2 )F, 



dz 



uv 



0. 



b.1 
2p 



The differential Eqs. (|5j) and (||) and the BC's (§D and 
wardly. The solutions are 



Fw(p,q) 



sin 



Z — ,u; 



x 



2sin(cu7r/2) VA' 



\p 
2q 

q 
p 



(58) 

can be solved straightfor- 

(59) 



1/2 



- p{-w)It 



uj/2 



±1 

2p 



where I± v are modified Bessel functions, and u is given by 

io = Vl -4A = v 7 ! - N a /N, N cr = 128/(3Cvr 2 ) 

Furthermore, 



and 



Z(q/A,u) 



71 



2 sin(^7r/2) 



\p(u)R (q/A, -co) - p{-u)R {q/A, uj)} 



(60) 

(61) 
(62) 



R{q/A,u) 



Iuj/2 



+ 4/2 



1 [q 
2Va 



A 



A X 1 . X 

-co S y---smy- 



A . /A 
2 Sm V2 



2a) 2 \l2A 1 ^ 2 



±1 
2 A 



(63) 
(64) 



The sinu;7r/2 results from the Wronskian between I^ (jl /2{x) and 1^/2 (x). 

By adopting the approximation (51) we have obtained an analytic expression for the 
Yukawa vertex IV Within this approximation, the a boson propagator Ag defined by 
Eq. (|40|) is related to T$ via Eq. (0). Such an expression is valid in the symmetric phase 
of the phase diagram. 



13 



IV. SCALING AND OTHER PROPERTIES 



In the previous section we have obtained nonperturbative solutions for the Yukawa vertex 
and scalar propagator within the ladder approximation. In this section we discuss some 
important properties of the Yukawa vertex and scalar propagator. 

Let us briefly state our objectives. First, we apply the Thouless criterion of the symmetry 
phase instability in order to derive the critical curve given in Eq. (||) . Subsequently, we show 
that near this curve the scalar propagator has a scaling form consistent with the general 
renormalization group theory of second order phase transitions. We find that the anomalous 
dimension of the propagator of the composite scalar fields is rj = 2 — uj with uj given by 
Eq. (pip. Moreover, we show that the Yukawa vertex has a scaling form consistent with 
power-law renormalizability. Second, we derive the peculiar behavior of lis near N = N cr . 
The phase transition at N = N CI is known as the CPT and is characterized by the absence 
of light unstable modes in the symmetric phase. Another characteristic feature of the CPT 
is the scaling law with essential singularity for the scalar boson and the fermion mass in the 
broken phase. This scaling law can be obtained by analytical continuation of lis in uj across 
the critical curve at N = N CT . 

In analogy with Ref. [^4j we investigate a few specific limits: 

(A) The large flavor limit (N — > oo), which means that the gauge interaction is negligible 
with respect to four-fermion interactions, i.e. A = 0, thus uj — 1. 

(B) Asymptotic or IR behavior of Ts(p + q,p) and As(q), i.e. p, g C A. 

(C) The behavior of lis at the critical coupling A = A c = 1/4, thus uj = 0. 

(D) The behavior of lis for A > A c , uj = iv, v = V4A — 1, i.e., analytic continuation across 
the critical curve at A = A c . 



A. Large flavor limit 



In the large flavor limit, the four-fermion interactions completely govern the dynamical 
breakdown of "chiral" symmetry. In this limit uj = 1 (A = 0), thus the Yukawa vertex (fH]) 
is Ts(p + q,p) = 1. Consequently, we obtain an expression for lis from Eq. ( p9|) by using 
Eq. ( |C13D and F\jy(p, q) = -FtrQo, q) = 1 at uj = 1. This leads to 



2NA 



71- 



1 



Aq q 2 
3A + 2A2 



(65) 



This expression is obtained by making use of the approximation (|C13|) . Naturally, the 
expression for 115(g) can be obtained by evaluating Eq. (^T|) with Y $ = 1- The result is 

J2 



n s (<z) 



2AA 



7i 2 q q- 
~8A + 3A2 



(66) 



see, e.g., Ref. and references therein [33]. Since only the first two terms on the right- 
hand side of Eqs. (|6S|) and (|66|) are important in the IR (q A), these equations differ 
about 10%. 
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B. Asymptotic behavior and scaling 



For values < uj < 1, the asymptotic behavior or IR behavior of T$ and Hg with 
(q/A) w ^> q/A can be derived by first considering the g<C A limit of Z, Eq. (|62j): 



7T 



2sin(cj7r/2) 



(i) 1/2 C Msi„h 



-In- + 5 u; 
2 g 



where 



6(u) 
C(u) 



1 p(^)(l + w)r(l + cj/2) w. A 
- In — — r — ; — in — . 

2 p(-u)(l -u)T(l -u/2) 4 8 



r(l + o;/2)r(l -u/2)' 
In this limit, the function Fjjy(p, q) with fermion momentum p = A can be expressed as 



(67) 

(68) 
(69) 



Fw(A,q) 



+ 



2u 



[1 - cothy) , 



1+u (l-uj 2 ) 
Thus, by using Eq. (Ek3), the asymptotic form for II5 reads 



V = -Jn- + d(wj 
2 q 



(70) 



n 5 (?) 



2NA 



7T^ 



+ 



1-cu 



2^2 



^1 — cothj/) 



Hence 

n s (?) 

where 



2NA 



7T 



A 



O {(q/A) 2 -) + {(q/A) 2 ) 



g< A, 



(71) 



(72) 



B(u) = 



16a; p(— u ) r (l — f ) 

;i + ^) 3 (i-cu) r(i + f) IV 8 



(73) 



One can show that -B(l) = 4/3, which is in agreement with Eq. (|B5|). The expression ( |72|) 
for the asymptotic behavior of Hs(q) is valid for < u < 1, but not for = (A = A c ). 
The inverse propagator A^ 1 that follows from Eqs. (|40D and ( [72]) is given by 



A* 1 (9) 



2B(u)NA 



7T- 



B(uj) \g g, 



(74) 



The instability of the symmetric phase is signalized by the vanishing of A^ x (g = 0). This is 
nothing else than the Thouless criterion for a phase transition of the second kind []33| which 
leads to the critical curve 



g = g c , < u < 1 (N > N CT ), g > -. 



(75) 
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Thus the curve g = g c is a line of UV stable fixed points. On the critical line the scalar 
propagator scales as 



A*(g) 



7T 



2B(lu)NA \ q 



A 



2-r? 



T] = 2-UJ, 



(76) 



where 77 is the anomalous dimension. 

On the other hand, one can see that on the line uj = (N = N cr ), g < 1/4, A^ x (g = 0) 
does not vanish. Nevertheless, as we shall show in Sec. [V|, this line is also the phase transition 
line but of a special type. 

The scaling form for Ts is obtained by considering only the leading term in Eq. 
Thus the Z function scales as 



Z(q/A,u) 



71 



2sin(cjvr/2) 2 T(l-u/2) \8 
In this way the Yukawa vertex can be written as 



-w/4 



q \ (i-w)/2 
A, 



(77) 



- A \ (ij-i)/a 

F s (p + q,p) ~ 1 ( - J [FiR(p/q)9(q - p) + F\jv(q/p)9(p - q)] 



(78) 



where, for p, q <C A, 



fa-l)/2 



Fm{p/q), F vw (p,q) « ( - 

9 



(»?-l)/2 



and 



Fm(p/q) 
Fw(q/p) 



2sin(c;7r/2) 2 r(l-w/2) / A 



7T p(o») (1+W) 

2 r(l-w/2) /A W4 



w/4 



sm 



\p 
2q 



(79) 



(80) 



p(cu) (l+a» 



1/2 



2p 



- p(-u)I u/2 



±1 

2p 



(81) 



An important consequence of the scaling behavior of the scalar propagator (Eq. (|76|)) and 
of the Yukawa vertex (Eq. (|78|)) is that, in using them, one finds that the four-fermion 
scattering amplitudes scale as 



Fs(pi + q,Pi)A s (q)T s (p2,P2 + q) oc -, pi, p 2 < q < A. 



(82) 



This scaling form reveals the long range nature and power-law renormalizability of the four- 



fermion interactions at the phase transition line 35 
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C. At the critical coupling 



At the critical value of A, i.e., uj = (A c = 1/4), we can derive in analogy with Ref. 
that for p ^> q 



(w = 0) F vv {p,q)^2 



A 



e 3 - 2 + ln(p/g) 
e 3 - ln(?/A) 



0{q 2 /p 2 Hq/p)) 



where 



1 



+ r [y/l/8) y/T/8 sin y/T/l 



e 2 = A ( ^/V 



^ U/V8 V^siny^ 



V / l78cos v /l78--sin v /lA 

/V8 cos - 2 sin 

e 3 = 2-7 + -ln2--, 

I ei 

with 7 the Euler gamma and K the modified Bessel function of the third kind. 
In the infrared, i.e., q <C A, Us can be written as 



{u = 0) n 5 (g) 



2AA 



7T 



4 
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+ 0(g 2 /A 2 ln(g/A)) 



ln(g/A) - e 3 

This straightforwardly follows from the insertion of Eq. fl83|) in Eq. 



(83) 

(84) 
(85) 
(86) 



(87) 



D. Analytic continuation across the critical curve 

Since the expression for the a boson vacuum polarization is symmetric under replacement 
of to by —uj, it can be analytically continued to the values A > A c . This holds in replacing 
uj by iv in Eq. ( pHj ) with Fuv given by Eq. (|5"0"D, where 

v = V4A - 1. (88) 

In the infrared (g <C A), it means that lis can be written as 

2/ = ^ln- + z/0(z/ 2 ), (89) 
z q 

where we have used Eq. ([H|) with uj — > iv, and where 0(z^ 2 ) = 5{iv)/iv. 

The four limits of rig described above are very useful for illustrating the resonance 
structure of the bound states and peculiar dynamics of the CPT, see Sec. |V|. 

To conclude this section let us mention that at zero a boson momentum (q = 0), we 
obtain 

r s (p,p) = F vv (p, q = 0) = — (|) . (90) 



2NA 



7P 



4(1 - v 1 



1 + f 



2\2 



1 + 1/ 



2^2 



cot y 
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V. LIGHT RESONANCES AND THE CONFORMAL PHASE TRANSITION 



In this section we analyze the behavior of the a boson propagator near the critical line 
in the symmetric phase (g < g c ), where the a and ir boson are degenerate. We will show 
that for u > (N > N cr ) the scalar composites (a and 7r bosons) are resonances (unstable 
modes) described by a complex pole in their respective propagators. The complex pole in As 
should lie on a second or higher Riemann sheet (i.e., not on the first (physical) sheet) of the 
complex plane of the Minkowskian momentum p 2 , because unitarity (causality) demands 
that A s{p) is analytic in the upper half of the complex p -pla.ne, where p^ is the "time" 
component of the Minkowski momentum p 2 = p\— p 2 . 

From Eq. (|74]) the complex pole can be computed. First we rotate back to Minkowski 



(91) 



space, p 2 — > p 2 M exp(— in). Subsequently, the complex poles are given by 



p 2 M = |m CT | 2 exp(-?6>), A/(p M ) 
The equation for the imaginary part reads 



2iVA 
n 2 g 



+ n s (p 



M) 



0. 







sin 



uie 



7T 



(92) 



with the solution 







-7T 



2nn 



(93) 



where n is an odd integer. Hence for values < lu < 1 it follows from Eq. ( p3|) that the 
complex pole does not lie on the physical sheet of p 2 .| Since cosu;(6' + 7r)/2 = —1, we find 
that the solution for \m„\ is 



A [g c gBiuj) 



1/uj 



^9 = 9c~ 9, 



(94) 



consequently the critical exponent v — 1/uj |^6|-0 Equation ( p4|) describes how the mass of 
the pole vanishes as g is tuned toward the critical line. 

The propagator As is of the form given by Eq. ([74]) and in Minkowski space, with the 
definition p = \fp 2 ~ , it can be written as follows: 



As(p) 



7T 



99c 



2NA Ag 



-1 



( _i)W2 ( p /\ mu \y 



(95) 



5 We denote the first (physical) Riemann sheet of p 2 by angles 9 with < — 9 < 2tt (the origin is 
a branch point with a branch cut along the positive real axis). 

6 The critical exponents v and rj coincide with those found in Ref. |p6| . Note, however, that in 
Ref. |3^] r\ was obtained assuming the validity of scaling relations between the critical exponents of 
the theory. Thus the independent computation of g in the present paper gives, in fact, a proof 
of the scaling relations. 
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with \m a \ given by Eq. (|94 

Re(A ff (p)) 
Im (A s (p)) 



Then, the real and imaginary part of As are 



7T 0& 



(p/\m a \) u cos <f\ 



2NAAg ) (p/\m a \) 2uJ + 2(p/\m r7 \) UJ cosy? + 1' 
f Jjf_99c\ (p/\m a \) u sirup 



\2NAAg J {p/\m a \) 2ul + 2{p/\m a \Y cos y? + 1' 



where = ixu)/2. The absolute value of the imaginary part has a maximum at p 
and the maximum is 



(96) 
(97) 
\m a \ 

(98) 

This shows that when g approaches g c from below (g | g c ), \m a \ goes to zero (|m CT | — > 0) 
and that the maximum of the absolute value of the imaginary part of As approaches infinity 
(-ImAsCKI) -> oo). 

We define a width over mass ratio r/lmJ as follows 



Im(A 5 (|m CT |)) 



7T 0& 



sin (f 



2NA Ag2[cos(p + l]' 



P+ 



m n 



P- 



Im(A s (p ± )) = ~Im(A s (|m CT |)). 



(99) 



Thus the width is the difference between the momenta at which — Im(As) equals 1/2 of the 
maximum value of — Im(As). Solving Eq. ( p9|) by making use of Eqs. ( |9*7|) and (p8|) gives 



2 + cos + a/(2 + cosv9) 2 - 1 



2 + cos if — a/(2 + cos<^) 2 — 1 



(100) 



Thus, as the mass scale of the pole is made small by approaching the critical line, the 
resonance is not described by a narrow Breit-Wigner type, because the width over mass 
ratio is rather large. Consequently, the resonance does not have the Lorentzian shape which 
is a characteristic feature of the Breit-Wigner resonance (note that even in pure NJL model 
(u = 1) the resonance is not narrow in contrast to four- dimensional NJL model). The 
above expression shows also that r/|m CT | increases when uo — > 0, and the resonance becomes 
broader. 

A description of the resonance structure is provided by a plot of lmAs(p). This is 
illustrated in Fig. |9| in which Im(A5(p))/As(0) is drawn as a function of the energy scale 
p/\m a \ for various values of to. 



A. Absence of light resonances near N cr 

The existence of light resonances whose mass vanishes as the transition is approached 
from the side of symmetric phase in (2+l)-dimensional theories is relevant for describing 



spin excitations in high-T c cuprate superconductors (see the paper by Kim and Lee [^5[ and 
references therein). Such resonances can be considered as precursors of the antiferromagnetic 
transition. It is known that QED3 by itself cannot give rise to light excitations in the 
symmetric phase p7|,|5|j6|. This is one of the main features of the so-called conformal phase 
transition: the absence of light excitations (composites) in the symmetric phase as the 
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transition is approached (in the broken phase massless "normal" Goldstone bosons appear). 
This unusual behavior can be attributed to the long range nature of the gauge interaction 
in the model under consideration. Another characteristic feature of the CPT is the scaling 
law with essential singularity for the dynamical fermion mass in the broken phase ||. 

From the side of the symmetric phase there is no sign indicating the occurrence of a 
phase transition. This means that the correlation length remains finite in the symmetric 
phase even close to the critical point (the Thouless criterion is not valid). In QED3 the CPT 
occurs at A = A c (N = N CT ) where the symmetry is dynamically broken by a "marginal" 
operator (a long range interaction). Though continuous, the CPT is not a second order 
phase transition. This is reflected by the singular behavior of some of the critical exponents 
(e.g., v and /3, see Ref. [[36|D as u goes to zero. The absence of a light complex pole in the 
a boson propagator illustrates the CPT in GNJL model in 2+1 dimensions. At u = the 
cr boson vacuum polarization is given by Eq. (|87|) in the infrared. If there has to be a light 
excitation in the symmetric phase then there must be a complex pole p 2 M = \m a \ 2 exp(— i9) 
in As with \m a \ <C A as g < 1/4. From Eq. fl37|), we then should find zeros of A^ 1 at 



\g J [ln(A/K|) + e 3 ] 2 + (^ + vr)V4' 
+ (102) 



For g < g c — 1/4, there are no solutions satisfying |m CT | <C A, hence if there is a pole it will 
be heavy, i.e., \m a \ ~ A. Therefore at A = A c and g < 1/4 there are no light resonances in 
the 2 + 1-dimensional GNJL model. 

What happens with the cx-boson propagator if we analytically continue it to the values 
X > X c , N < N cr 7 By doing so, we remain in the massless chiral symmetric phase, but we 
just end up in the "wrong vacuum" (the chiral symmetric solution becomes unstable).]] The 
7r and a bosons are tachyons for such a solution. Thus the border of the stable symmetric 
solution A = A c , g < 1/4 is also the phase transition line. 

Let us show that there are indeed tachyons (with imaginary mass m 2 < 0) when A > A c . 
For this we need to show that A5 has a real pole in Euclidean momentum space. Assuming 
that the pole lies in the infrared, \m a \ <C A, we can use Eq. (f^TJ), where u has been replaced 
by iv, v given by Eq. fl88|). The tachyonic pole is given by 



2iVA 
n 2 g 



ll s (\m a \)=0. (103) 



From this we derive the solution 



\ m cr\ ( 277, 7T 2(3 t 2 



A V v v 

where 



exp ^ + 20(z/) , (104) 



7 A phase transition is by definition described by a nonanalytic behavior in the theory parameters 
(coupling constants, temperature, etc.). Therefore by analytical continuation one cannot go from 
one phase into another. 
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p = tan" 1 ^ -. (105) 

g-2X(g + \) 1 ; 

As is well established, the tachyon with the largest \m a \ in the physical region corresponds 
to n = 1. 

It is clear that the tachyons in the symmetric solution appear also when we cross the 
upper part (A < A c , g > 1/4) of the critical curve. However, the difference between this part 
of the critical line and the line A = A c , g < 1/4 is that we have light composites (resonances) 
near the first line while they are absent near the last one. 

If we now consider the limit A [ A c , i.e., we approach the transition from the side of the 
broken phase, we obtain the scaling law with essential singularity, 

" exp f - + 2(f)(0) j exp f ^ ) . (106) 



A " v V V v/4A 

This scaling law with essential singularity is obtained by analytical continuation of the 
solution in the symmetric phase (A < A c ) to the broken phase (A > A c ). Thus the tachyonic 
(unphysical) solution in the broken phase leads to a scaling law which is proportional to the 
scaling law given by the fermion mass and a boson mass in the broken phase M. 



VI. CONCLUSION 

In this paper we studied a nonlinear equation for the running coupling in QED3 which 
can be considered as the analogue of the ladder approximation for the fermion propagator. 
We solved our equation both analytically and numerically. We found that the vacuum 
polarization operator, obtained through a nonperturbative solution of the equation, has the 
same infrared asymptotics as the one-loop expression: Tl(p) ~ Ca/p, C ~ 1 + 1/14AT. Thus, 
we have showed that a nontrivial IR fixed point persists in the nonperturbative solution. 
Moreover, quantitatively there is only slight difference between the one-loop result and the 
nonperturbative solution even at the number of fermions N = 2. 

We then proceeded with studying the GNJL model in the symmetric phase with massless 
fermions. We solved an equation for the Yukawa vertex in the approximation where the 
full Bethe-Salpeter kernel is replaced by the planar one photon exchange graph with bare 
fermion-photon vertices. The obtained Yukawa vertex was used for calculation of the scalar 
composites propagator. The phase transition curve was determined from the condition 
of instability of the symmetric solution. We established the existence of light excitations 
(resonances) in the symmetric phase for values of A^ > N cx ^ 4 (A < A c ), provided the 
four- fermion coupling (g > 1/4) is near its critical value along the critical curve (§). As 
g < 1/4 and A" approaches Ac r from above the light excitations are absent and the situation 
resembles pure QED 3 . 

The field theoretical models, like QED 3 and the GNJL model, often appear in the long- 
wavelength limit of microscopic lattice models used for description of high-T c samples. For 
instance, in a spin-charge separation Ansatz for the t — J model, where spin is described by 
fermionic spinons and charge is described by bosonic holons (or vice versa), a "statistical" 
U(l) gauge interaction appears naturally in the theory along with four-fermion interactions 



see, for example, Refs. [ll], [20], and f2"B|). It was argued in Ref. [25] that QED 3 , with 
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fermions treated as spinons, might serve as a possible candidate for describing the undoped 
and underdoped cuprates. For physical N(= 2) the chiral symmetry broken phase of QED 3 
(with a dynamical mass generation) should correspond to an antiferromagnetic ordering in 
undoped cuprates pSj , while the symmetric phase (for larger N) would describe some kind 
of a spin liquid. 

Recently spin excitations (particle-hole bound states) have been observed in the normal 
state (and in the superconducting state) of underdoped and optimally doped cuprates such 
as YBa 2 Cu30 6+:r and La 2 _ a; Sr 3 .Cu04, where x is the amount of doping, see Ref. |M| and 
references therein. The dynamic susceptibility x" describing antiferromagnetic correlations 
near wave vector Q = (it, tc) has a broad peak whose energy comes down as the doping is 
reduced. The height of the peak increases as the doping is reduced and the antiferromagnetic 
transition approached. 

As was proposed in Ref. J25|, QED 3 could describe these particle-hole excitations. How- 
ever, from the point of view of the present paper, pure QED 3 cannot be applied for describing 
such spin excitations because of absence of light resonances in the symmetric phase of the 
model. In our opinion, the GNJL model serves this purpose better since light excitations 
appear near the critical curve (|3|) on both sides. Moreover, the mass of resonances decreases 
as the phase transition is approached (along the trajectory N, or A, is fixed and g | g c ) while 
their peaks become sharper as g approaches g c . All these features are in qualitative agree- 
ment with the experimental picture if we assume that the four-fermion coupling g depends 
on the doping in such a way that g increases when the doping is reduced. 

A problem, however, is that, in case of cuprate superconductors, the physically relevant 
number of flavors equals two (N = 2) which is less than iV cr ~ 4. This means that one 
would get the dynamically broken symmetry phase corresponding to the Neel ordered state 
at any doping in both QED 3 and the GNJL model. Kim and Lee |2^| proposed a mechanism 
to lower iV CI . (and make N cr < 2) in pure QED 3 by taking into account the effect of doping 
which screens the time-component of the gauge field and halves N CT , due to additional 
coupling of the gauge field to charged scalar fields. However, another way out of such 



a dilemma appears if we invoke the arguments of Appelquist et al. [40| that ladder SD 
equations usually overestimate the critical value N cr . These authors suggest that the true 
critical value is iV cr = 3/2. Thus for the physical case of N = 2 the spontaneous symmetry 
breaking does not occur and the system will be in the symmetric phase when the doping 
exceeds some critical value. It would be quite interesting to find out a truncated set of SD 
equations giving such a small critical N CT . 
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APPENDIX A: BOX DIAGRAM 



In this appendix we compute the box function B of Eq. (PJ). We start by contracting 
the 7's in Eq. ([15]) and evaluate the traces. The result is 



B(p\k\p-k) = i 



d 3 r 



M 



(2tt) s 



[bi(p, k , r ) + hip, k, r) + hip, k, r)] , 



(Al) 



where 



hip, k,r) 



— 16(A; • r) 2 — Ak ■ r(4fc • p — p 2 ) + 4(/c 2 — Ak ■ p)p ■ r — 24k ■ rp ■ r — 16(p • r) 



-24k ■ rr 2 + 4(A; 2 - 3k ■ p + p 2 )r 2 - 24p ■ rr 2 - 12r 4 
1 



x 



(r + p) 2 (r + p + /c) 2 (r + fc) 2 r 2 



(A2) 



and 



hip, k,r) 



Ak ■ rp 2 + 4(2A; • p + p 2 )p ■ r + 8(p • r) 2 + 4A; • rr 2 + 4(2A; • p + p 2 )r 2 

1 



hip,k,r) 



+ 12p-rr z + 4r q , 

r A \r + p) 4 (r + p + k) A 

8k ■ rp ■ r — 4k ■ pr 2 + 4k ■ rr 2 + 4p ■ rr 2 + 4r 4 

r 4 (r + p) 2 (r + A;) 2 



(A3) 
(A4) 



The traces have been performed with the help of FeynCalc jO]. Subsequently, we cancel the 
r dependence in the numerators of Eqs. ( [A2| )-(|A"4]) without shifting the integration variable. 
In this way the box function B can expressed as follows (in Euclidean formulation): 



Bi-p 2 ,-k 2 ,-p-k) 



d 3 7 



(2*) 



+ 



1 



+ 



1 



+ 



2A; 2 



r 2 (r + k) 2 r 2 ir + k + p) 2 (r + p) 2 (r + A;) 2 

k 2 4k 2 
+ 



V+p) 4 r 2 (r + p) 2 (r + k) 2 (r + p) 4 (r + k + p) 2 (r + + p) 2 (r + p) 2 (r + A;) 2 
2A; 4 2fc-p 2fc-p 



+ 



r 2 (r + k + p) 2 (r + p) 2 (r + A;) 2 r 2 (r + p) 2 (r + A;) 2 r 2 (r + p) 2 (r + A; + p) 2 
2k ■ p 4k ■ p 4k 2 k ■ p 



r 2 (r + k + p) 2 ir + A;) 2 (r + k + p) 2 (r + p) 2 (r + A;) 2 r 2 (r + A; + p) 2 (r + p) 2 (r + k)' 



+ 
+ 



2A; • r 



4/c • r 



+ 



P 



P 



P 



r 4 (r + A;) 2 r 2 (r + p) 2 (r + A;) 2 r 4 (r + A;) 2 r 2 (r+p) 4 r 2 (r + A; + p) 2 (r + Ac) 2 
4p2 ^,2^2 5/c 2 p 2 



+ 



(r + A: + p) 2 ir + p) 2 (r + k) 2 r 2 (r + p) 4 (r + k + p) 2 r 2 (r + A; + p) 2 (r + p) 2 (r + k)' 



4k ■ pp 2 



2k ■ rp 2 



p 



r 2 (r + k + p) 2 (r + p) 2 (r + A;) 2 r 4 (r + p) 2 (r + A;) 2 r 4 (r + p) 2 (r + A;) 2 
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2p 4 



4p ■ r 



2p ■ r 



r 2 (r + k + p) 2 (r + p) 2 (r + k) 2 r 2 (r + p) 2 (r + k) 2 r 2 (r + k + p) 2 (r + k) 2 

2p 2 p ■ r 2r 2 
r 4 (r + p) 2 (r + k) 2 (r + k + p) 2 (r + p) 2 (r + A;) 2 



(A5) 



For the explicit calculation of the integral, only a handful of integrals are involved. These 
integrals are 



d 3 r 
(2^)3 



k 2 



d 3 r 



2^2 



k 2 p 



r 4 r 4 (k + r) 2 
k 2 



0. 



E (2tt) 3 [r 4 (k + r) 2 (p + r) 2 r 4 (k + rf 

d 3 r 1 



fc • p 



, E (2tt) 3 r 2 (A; + r) 2 
d?r 1 



1 

8k' 

1 



(27r) 3 r 2 (r + A;) 2 (r + p) 2 8kp\k — p| ' 



ci 3 ? 



(27r) 3 r 2 (r + fc) 2 (r + p) 2 (r + k + p) 2 8k ■ p 



kp\k — p\ kp\k + p\ 



(A6) 
(A7) 
(A8) 
(A9) 
(A10) 



where \k — p\ = \/Jk — p) 2 . In the computation of the last integral ( A10|) , we first have 
rewritten the left-hand side as the sum of four so-called triangle diagrams (diagrams of the 
form of Eq. (JA9|)), the subsequent integration is then straightforward. 
The final result reads 



B(-p 2 ,-k 2 ,-p-k) = - + __ 



k ■ p 



+ 



k ■ p 



k p Akp\k + p\ Akp\k — p 



(2k 2 + kp + 2p 2 
2kp\k + p\ 



(2k 2 + kp + 2p 2 ) (2k 4 + 5k 2 p 2 + 2p 4 ) + (2A; 4 + 5k 2 p 2 + 2p 4 ) 



2kp\k — p\ 



A(k ■ p)kp\k + p\ 



A(k ■ p)kp\k — p\ 



(All) 



If an additional integration over the angle between p and k follows we can simplify this 
equation because of the symmetry p ■ k — > — p ■ k. However, if we make use of this symmetry, 
we should take the principal value for angular integration, since the 1/p ■ k singularity no 
longer explicitly cancels. 

At various places, we need the following angular integrals:^] 



1 



1 



An \k — p\ m&x^jp)' 

dQ 1 1 , k+p 

m 



dQ k ■ p 



k 2 p 2 



4ir \k — p\ 3 max(/c 3 ,p 3 



47r (k + p) 2 2kp \k — p 



(A12) 
(A13) 



? The angular measure is f dQ = 2ir d6 sin 6, with cos 9 = k • p/kp. 
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dQ k ■ p 1 



An (k + p) 2 2 Akp 
dVt (k-pf _(k 2 + p 2 ) 
4^ (k + p) 2 " 4 

and the Cauchy principal value integral 

1 1 



(k 2 +p 2 ) j k + p 



\k — p 
(k 2 +p 



2kp 



2 V k + p 
In 



(A14) 
(A15) 



In 



+ p + \/k 2 +p 2 



4tt k-p\k-p\ kp^Jk 2 + p 2 \\k 



P\ 



Vk^+ 



P 



(A16) 



Note that we use the notation k = yk 2 , p = sjp 2 for the scalar quantities in the right-hand 
side expressions of Eqs. (|A12|) -( [A1~6| ). 

Here, we compute also a general form of the integral over the kernel K(p, k) given by 
Eq. (|D, i.e., 



1(6) 



dk —K(p, k) 



p" 



dtr 



lit 

~6~ 



1 2t 4 + 5t 2 + 2 . 
H : — , smh t 



+ 



i 

p/A 
1 



dtt 



-is 



t 2t 2 VT+¥ 

lit 1 2t 4 + 5t 2 + 2 



6 



o 



dt (t' + t- 1 - 6 ) 



t 2t 2 VT+t 2 

lit 1 2t 4 + 5t 2 + 2 

2t 2 VT+¥ 



sinh 1 1 

sinh -1 1 



2p„ 
H hm t 

15 A t^o 



2-5 



+ o 



p 
A 2 



(A17) 



with —3 < 5 < 2. Exact solutions exist when 6 is an integer. In that case, one can make 
the transformation t —>■ sinh Ins = (s 2 — l)/2s, after which the integral can be written as a 
sum of Spence functions. The result is 



1(0) 



7T 
1 



2 5 

2' 



7T 



23 
18' 



1(2) 



64 ~ 4' 



(A18) 



APPENDIX B: DERIVATION OF THE NONLOCAL GAUGE IN THE GNJL 

MODEL 

In this appendix we derive the nonlocal gauge £ in the GNJL model in order to set A = 1, 
so that the WTI is satisfied if one uses the bare vertex approximation. For generality we 
consider the case of arbitrary dimensions D and in presence of the mass function B. 

We introduce the nonlocal gauge function ^(k 2 ) by writing the full photon propagator in 
the form 

e 2 D, v (k) = - [g, v - C(* 2 )^) (Bl) 
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where d(k 2 ) = e 2 /(l + I1(A; 2 )), ((k 2 ) = 1 — £,{k 2 ). The SD equation for the fermion wave 
function renormalization A is given by 



A(p 2 



1 + 



d D q A(q 2 



p 2 J (2tt) d q 2 A 2 (q 2 ) + B 2 (q 2 ) 



X 



d(k 2 



(D — 2)p ■ q + ip-q 



2(p 2 q 2 - (p • q) 2 ) 



k 2 



C(k 2 



-p-q [A s (k 2 ) + A P (A; 2 )]|, k=p-q. 



^2) 



Introducing the variables k 2 = x + y — 2^/xy cos#, x = p 2 , y = q 2 , and performing the 
integration over all angles except the angle 8, we get 



A 2 



p 2 (A(p 2 ) - 1) = C D J dy 



y {D - 2)/2 A{y) 



yA 2 (y) + B 2 (y) 
o o 



d6 sin D - 2 e 



x <^ d{k z ) 



xy cos 6(D -2 + ((k 2 )) sin 2 6 ^ 
p 2xy—((k ) 



xycosO [A 5 (A; 2 ) + A P (k' 



(B3) 



where C^ 1 = 2 D n ( - D+1 ^ 2 T((D — l)/2). Following the works by Kugo et al. and Kondo et al. 
|pE2]| (see also Ref. ||43|| ) we perform now the 9 integration by parts in terms containing the 
first power of cos#: 



p 2 (A(p 2 ) - 1) 



x 



D — 1 



yD-l 



yA 2 (y) + B 2 (y) 



(z D - 2 d(z)C(z))' — (D — 2)z D ~ z (d(z) - zd\z)) 



[A s (z) + A P (z)}' 



(B4) 



where the prime denotes the differentiation with respect to z = k 2 . 

The requirement A(p 2 ) = 1 is fulfilled by choosing ((z) such that the expression in curly 
brackets in Eq. (|B4]) vanishes. This gives the first order differential equation for Q{z) which 
is easy to integrate, 



D-2 

z D - 2 d{z) 



dtt D ~ 3 [d{t)-td'{t)\ + 



z D - 2 d(z) 



dtt D - 1 [A s {z) + A P {z)}' (B5) 



: ) I „ = in order to eliminate 

2=0 



(the integration constant was fixed by requiring \z D 2 d(z)£( 

the singularity at z = in Ci z ))- The last equation finally leads to the following expression 
for f (z): 
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M= D - i - {D :£%) 2) i* tD ~' m 



z 

^- ) jdtt D - 1 [A s (z) + A P (z)]'. (B6) 



z D-2, 



For D = 3 we take 



e 2 



d ^ = T+WSj ~ NC k ' k<<a - (B7) 

with k = y/k?, and assume the following form for scalar propagators in the symmetric phase 
and near the critical line (see Eq. (|76|)) 

a /A 2 V 

A s (k) = A P (k) = -- - , (B8) 



A \k\ 

where a is some constant and the power < 7 < 1 (Eq. ( p8|) is verified a posteriori when 
solving the SD equation (|40| ) for the scalar propagator). We obtain, from Eq. 



m *^rpr\ (B9 ) 



3 4(2 -7) V k 



In absence of the four-fermion interaction we get the famous nonlocal gauge £ = 2/3 [|T6 



Carena et al. p]J have included exchanges by the bare scalar propagators what corresponds 
to taking 7 = 1/2, a = 4/N {C = 1 in their leading order of the 1/N approximation for the 
photon vacuum polarization). Equation ( |B9| ) then gives £ = 1/3 in accordance with their 
findings. 

Our Eq. (74) for the scalar propagator gives the exponent 7 = u/2 and contribution due 
to the exchange of scalars into becomes suppressed (since to < 1) and we are left with 
Nash's nonlocal gauge £ = 2/3. 



APPENDIX C: APPROXIMATION FOR THE YUKAWA VERTEX 

As was mentioned in Sec. |T|, in order to resolve the angular dependence of the Yukawa 
vertex function Fx, we expand it, together with the kernels of Eqs. (fl8D and (|49|), in Legendre 
polynomials P n . We write 

00 

= ^Y n (fc,p)P n (cos/3), (CI) 

n=0 

(C2) 

where cos a = p ■ q/pq, cos/5 = p ■ k/pk, and cos 7 = q ■ k/qk. The Legendre polynomials P n 
satisfy 



Fx(p + q,p) = ^2fn(p,q)P n (cosa), 

n=0 

(fc + g )2 = 22an{k,q)P n (cosj), 



n=0 
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1 

— P m (cosa)P„(cosa) = - / dx P m (x)P n (x) = mn . (C3) 
47T Z J In + 1 



-i 



With the above defined expansions, and by making use of the identity 

1 f n /*27T £ 

— / da sin a / d9 P n (cos a) Pi(cos (3) = — — ^ — p(cos7), (C4) 
47T J Jo 21 + 1 

where cos (3 = cos a cos 7+ sin a sin 7 cos 9, Eq. ([48]) for the Yukawa vertex can be represented 
as the set of equations for harmonics f): 

/'^ 00 00 

dkNi(k,p) Y] ^ Q mn a m (k, q)f n (k, q), (C5) 

m=0 n=0 



where 

= 5 / dxP,(x)P m (x)P ri (x) = n /' ) ^ ( ^- m( ^T" g ' , (C6) 

2,/ (s - l)\ [s - my. [s - n)\ y~) s (2s + 1) 

where 2s = I + m + n and (a)& = T(a + k)/T(a). The coefficients C/ mn are zero unless 
I + m + n = 2s is even and a triangle with sides /, m, n exists, i.e., \l — m\ < n < I + m.f\ 
Furthermore, Eq. ( (4"5| ) can be written as 



A 

00 



2N f j, a n(k,q)f n (k,q) 



^ ^ — n 







2n + l 

n=0 



(C7) 



Within the approximation (51), i.e., keeping the zero-order term f only in right-hand 



sides of Eqs. (|C!5D and (|G7|) , we get 



A 



/o(p,g) = ! + dkN (k,p)a (k,q)f (k,q), (C8) 


A 

2iV /" 

n 5(?) = ^/ dka (k,q)f (k,q). (C9) 



The functions iVo and ao are straightforwardly obtained from inverting Eqs. (|C1| ) and flC2]). 
This gives 

Wo(M = !(^ + !<ez*>, (cio) 



9 We thank L. P. Kok for pointing out the paper by Askey et al. [31]. 



28 



and 



/dtt (k 2 + q- k) 
4^ (fc _|_ g )2 = aiR ^' - fc ) + a uv(fc, - <?), (Cll) 

/, n 1 (k 2 -q 2 )k + q 1 (/c 2 - g 2 ) fc + g 



In order to be able to solve the equations for F IR and Fuv given by Eq. Q5^), we approximate 
the functions aiR and auv as follows: 



& 2 q 2 1 

9) ~ , «uv(fc, g) ~ 1 - ^2 , %(?, 9) = auv(?, g) = ( C13 ) 



The validity of this approximation is addressed in Sec. [TV A . 

The lowest order harmonic /o of the particular Legendre expansion given in Eq. 
is expressed in terms of the so-called infrared (IR) function F IR and the ultraviolet (UV) 
function Fuv, see Eq. Q52]). These functions describe the following asymptotic behavior of 
the Yukawa vertex: 

\imT s (p + q,p) = llimFuv(p,g), (C14) 
UmT s (p + q,p) = 1 lim Fni(p, <?)• ( C15 ) 

<?>P <?>p 

The fact that both these asymptotic limits of Ts are described by fo through Fir and 
F\jy guarantees the validity of the approximation flSip. This crucial point is explained 
in more detail in Ref. [^4j], where this approximation is referred to as "the two-channel 
approximation." 

Then, by making use of Eqs. fl52|), (|C10|) , and ( |C13| ) we get 

V q 

/k f k 

dk—F m (k } q) + \ J dk—F m (k } q) 

p 

A 

+x J dk U i -&) fmk g) ' (ci6) 



q 

g 



(p > q) F uv (p, q) = 1 + A J dk £^F 1R (k, q ) + xjdkUl-^\ F uv (fc, q) 





A 



+ A / dfci^l-ijj?uv(*,g), (C17) 

and for the scalar vacuum polarization ( |C9| ) we can derive Eq. (0). The integral Eqs. ( p!6|) 
and ( pi7|) are equivalent to the second order differential equations given in Eqs. (0) and 
(j55|), with the four boundary conditions (|56p and 
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FIGURES 
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FIG. 1. Exact SD equations for the gauge boson propagator D^ u , the fermion propagator S, 
and the vertex T M . 




FIG. 2. The l/N expansion of the two-fermion, one-photon irreducible fermion-fermion scat- 
tering kernel. 




FIG. 3. SD equation for fermion propagator and vertex up to order l/N 2 . 
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FIG. 4. Closed SD equation for the gauge boson propagator in next-to- leading 1/N expansion. 
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FIG. 5. The box diagram B^ ua . 
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FIG. 6. Numerical solutions of Eq. (26) 
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FIG. 7. The SDE for the scalar propagator A s (p). 




FIG. 8. The SDE for the scalar Yukawa vertex T$ in the ladder approximation. 
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FIG. 9. The response function Im(A s (p))/A s (0) vs. p/\m a \ for u = 1.0, u = 0.8, u = 0.6. 



34 



